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Q-i Abstract 

We evaluate two-body decay modes of neutralinos in the Minimal Supersymmetric 
Standard Model with complex parameters (cMSSM). Assuming heavy scalar quarks 
we take into account all two-body decay channels involving charginos, neutralinos, 
(scalar) leptons, Higgs bosons and Standard Model gauge bosons. The evaluation of 
the decay widths is based on a full one-loop calculation including hard and soft QED 
radiation. Of particular phenomenological interest are decays involving the Lightest 
Supersymmetric Particle (LSP), i.e. the lightest neutralino, or a neutral or charged 
Higgs boson. For the chargino/neutralino sector we employ two different renormaliza- 
tion schemes, which differ in the treatment of the complex phases. In the numerical 
analysis we concentrate on the decay of the heaviest neutralino and show the results 
in the two different schemes. The higher-order corrections of the heaviest neutralino 
decay widths involving the LSP can easily reach a level of about 10—15%, while the cor- 
rections to the decays to Higgs bosons are up to 20—30%, translating into corrections of 
similar size in the respective branching ratios. The difference between the two schemes, 
indicating the size of unknown two-loop corrections, is less than 0(0.1%). These cor- 
rections are important for the correct interpretation of LSP and Higgs production at 
the LHC and at a future linear e + e~ collider. The results will be implemented into 
the Fortran code FeynHiggs. 
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1 Introduction 



One of the important tasks at the LHC is to search for physics beyond the Standard Model 
(SM), where the Minimal Supersymmetric Standard Model (MSSM) [lj is one of the leading 
candidates. 

Two related important tasks are investigating the mechanism of electroweak symmetry 
breaking, as well as the production and measurement of the properties of Cold Dark Matter 
(CDM). The Higgs searches currently ongoing at the LHC (and previously carried out at the 
Tevatron [2] and LEP [3]) address both those goals. The spectacular discovery of a Higgs-like 
particle with a mass around M# ~ 125 GeV, which has just been announced by ATLAS 
and CMS [lj, marks a milestone of an effort that has been ongoing for almost half a century 
and opens a new era of particle physics. Both ATLAS and CMS reported a clear excess 
around ~ 125 GeV in the two photon channel, as well as in the ZZ^ channel, whereas 
the analyses in other channels have a lower mass resolution and at present are largely less 
mature. The combined sensitivity in each of the experiments reaches about ~ 5er. The 
discovery is consistent with the predictions for the Higgs boson in the SM [5j, as well as with 
the predictions for the lightest Higgs boson in the MSSM j5j[6]. The latter model also offers 
a natural candidate for CDM, the Lightest Supersymmetric Particle (LSP), i.e. the lightest 
neutralino, Xi 0- Supersymmetry (SUSY) predicts two scalar partners for all SM fermions 
as well as fermionic partners to all SM bosons. Contrary to the case of the SM, in the MSSM 
two Higgs doublets are required. This results in five physical Higgs bosons instead of the 
single Higgs boson in the SM. These are the light and heavy CP-even Higgs bosons, h and 
H, the CP-odd Higgs boson, A, and the charged Higgs bosons, H ± . In the MSSM with 
complex parameters (cMSSM) the three neutral Higgs bosons mix [8HTT]. giving rise to the 
CP-mixed states hi, hi, h%. 

If SUSY is realized in nature and the scalar quarks and/or the gluino are in the kine- 
matic reach of the LHC, it is expected that these strongly interacting particles are copiously 
produced. The primarily produced strongly interacting particles subsequently decay via cas- 
cades to SM particles and (if i?-parity conservation is assumed, as we do) the LSP. One step 
in these decay chains is often the decay of a neutralino, Xi 3 4 , to a SM particle and the LSP, 
or as a competing process the neutralino decay to another SUSY particle accompanied by a 
SM particle. Also neutral and charged Higgs bosons are expected to be produced this way. 
Via these decays some characteristics of the LSP and/or Higgs bosons can be measured, see, 
e.g., Refs. [12j[T3] and references therein. At any future e + e~ collider (such as ILC or CLIC) 
a precision determination of the properties of the observed particles is expected [Lfp5]. (For 
combined LHC/ILC analyses and further prospects see Ref. [16].) Thus, if kinematically ac- 
cessible, the pair production of neutralinos with a subsequent decay to the LSP and/or Higgs 
bosons can yield important information about the lightest neutralino and the Higgs sector 
of the model. 

In order to yield a sufficient accuracy, at least one-loop corrections to the various neu- 
tralino decay modes have to be considered. In this paper we evaluate full one-loop corrections 
to neutralino decays in the cMSSM. If scalar quarks are sufficiently heavy (as in many GUT 
based models such as CMSSM, GMSB or AMSB, see for instance Ref. [TTJ ) a neutralino 
decay to a quark and a scalar quark is kinematically forbidden. Assuming heavy squarks we 
calculate the full one-loop correction to all two body decay modes (which are non-zero at 
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the tree-level), 



(2 = 2,3,4; j = 1,2) , 
(z = 2,3,4; j = 1,2) , 
(2 = 2,3,4; j<z; fc = 1,2,3) , 
i = 2,3,4; j < i) , 
{i = 2,3,4; £ = e,/i,r; fc = 1,2) , 



(1) 
(2) 
(3) 
(4) 
(5) 
(6) 



(2 = 2,3,4; £ = e,/i,r) . 



The total width is defined as the sum of the channels flTJ to (jB]), where we neglect the decays 
to colored particles as these will not be kinematically allowed for the scenarios considered in 
the numerical analysis. It should be noted that several modes are closed for nearly the whole 
MSSM parameter space due to the structure of the chargino and neutralino mass matrices 
(see below). Therefore, while we have evaluated analytically all neutralino decays, in our 
numerical analysis we will concentrate on the decays of the heaviest neutralino, xl- 

As explained above, we are especially interested in the branching ratios (BR) of the 
decays involving a Higgs boson, Eqs. (jTJ), (j3J) as part of an evaluation of a Higgs production 
cross section, and/or involving the LSP, Eqs. (131), (j4"|) as part of the measurement of CDM 
properties at the LHC, the ILC or CLIC. Consequently, it is not necessary to investigate 
three- or four-body decay modes. These only play a significant role once the two-body modes 
are kinematically forbidden, and thus the relevant BR's are zero. The same applies to two- 
body decay modes that exist only at the one-loop level, such as \°i ~~ >* X?7 ( see > f° r instance, 
Ref. [18J). While this channel is of 0(a 2 ), the size of the one-loop corrections to Eqs. (p"} to 
([6]) is of 0(a). We have numerically verified that the contribution of T(x° — y x?t) to the 
total width is completely negligible. 

Tree-level results for the neutralinos decays in the MSSM were presented in Refs. *T8Tf20] . 
The code SDECAY [21] includes all two-body decays of neutralinos at tree level. Tree-level 
studies of neutralino decays have shown that they could be invaluable in distinguishing 
between different patterns of supersymmetry breaking [22], as well as in detecting CP vio- 
lating effects at a linear collider [23TT28] or a muon collider [29]. Higher-order corrections to 
neutralino decays have been evaluated in various analyses over the last decade in the real 
MSSM (rMSSM), for which the on-shell renormalization of the chargino- neutralino sector 
was developed in Refs. [30-39J. In Ref. [40J three-body decays into the LSP and quarks 
were calculated including corrections to the masses of third generation fermions and SUSY 
particles. In Ref. [36], decays of the next-to-lightest neutralino to the lightest neutralino and 
two leptons were calculated at one-loop. The one-loop electroweak corrections to all two- 
body decay channels of neutralinos, evaluated in an on-shell renormalization scheme, have 
been implemented in the code SloopS [41] . Radiative corrections to a number of neutralino 
decay channels were also recently studied in Ref. [42] for the case of real parameters. A full 
one-loop calculation of the electroweak corrections to the partial width of the decay of a neu- 
tralino into a chargino and a W boson in the MSSM and NMSSM is presented in Ref. [43J, 
and made available with the code CNNDecays. In the cMSSM, the on-shell renormalization 
of the chargino-neutralino sector was first studied in Ref. [37J, and subsequently in Ref. [44J, 
and decays of type r(%° — > xjH^) were studied in Refs. [371138]. The approach to the renor- 
malization of the complex parameters differs slightly from that used in Ref. [45114*6] . where 
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chargino decays in the cMSSM at the one-loop level were analyzed. One important part of 
this work consists of a comparison of these two schemes. 

In this paper we present for the first time a full one-loop calculation for all non-hadronic 
two-body decay channels of a neutralino, taking into account soft and hard QED radiation, 
simultaneously and consistently evaluated in the cMSSM. The calculation is based on two 
independent set-ups that differ slightly in the inclusion of higher-order corrections to quan- 
tities used in one-loop corrections, i.e. in effects beyond the one-loop level. The two set-ups 
furthermore employ renormalization schemes in the chargino/neutralino sector that differ in 
their treatment of complex phases [37J[38j|l5j|l6] . The numerical results are shown for both 
set-ups and the (small) differences indicate the size of theoretical uncertainties beyond the 
one-loop level. 

The paper is organized as follows. In Sec. [2] we review the relevant sectors of the 
cMSSM and give all the details about the two different renormalization schemes in the 
chargino/neutralino sector. Details about the calculation can be found in Sec. [31 The nu- 
merical results for all decay channels are presented in Sec. HI The conclusions can be found 
in Sec. [5j The evaluation of the branching ratios of the neutralinos will be implemented into 
the Fortran code FeynHiggs [4TH5D] . 

2 The relevant sectors of the complex MSSM 

All channels (PQ) - (jH]) are calculated at the one-loop level, including real QED radiation. This 
requires the simultaneous renormalization of several sectors of the cMSSM. In the following 
subsections we introduce our notation for these sectors. Details about the two renormaliza- 
tion schemes used in the chargino/neutralino sector are given. The renormalization of the 
other sectors can be found in Refs. [45 1 146 1 151] . 

2.1 The chargino/neutralino sector of the cMSSM 

While many details about the renormalization of the cMSSM can already be found in 
Refs. [33I3B1S51B6J , we repeat here the most relevant aspects in order to give a complete 
picture and to facilitate the comparison between the two employed renormalization schemes. 
The chargino/neutralino sector contains two soft SUSY-breaking gaugino mass parameters 
Mi and M 2 corresponding to the bino and the wino fields, respectively, as well as the Higgs 
superfield mixing parameter fi, which, in general, can be complex. Since not all the possible 
phases of the cMSSM Lagrangian are physical, it is possible (without loss of generality) to 
choose some parameters real. This applies in particular to one out of the three parameters 
Mi, M2, and M3, the gluino mass parameter. For the numerical analysis in Sec. |4]we choose 
M2 to be real, however for the renormalization scheme I introduced below [45,46j, we do not 
make such an assumption, and the analytical derivation of the renormalization constants is 
performed for a complex M 2 , discussed further in Sec. HI 

The starting point for the renormalization procedure of the chargino/neutralino sector is 
the part of the Fourier transformed MSSM Lagrangian which is bilinear in the chargino and 
neutralino fields, 

4- x° = ^-xl + Tii u+xi - W [v*x T u^ u -xj - x7 [ux* V T ] 2J u +X j 
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+ 2 (xltfu-xl,+Xktfu+X°k 



Xl [N*YN t ] 



x£|NY*N T ]«a; + x?) , (7) 



already expressed in terms of the chargino and neutralino mass eigenstates xT an d Xki 
respectively, and i,j = 1,2 and fc, I = 1,2,3,4. The mass eigenstates can be determined 
via unitary transformations where the corresponding matrices diagonalize the chargino and 
neutralino mass matrix, X and Y, respectively. 

In the chargino case, two 2x2 matrices U and V are necessary for the diagonalization 
of the chargino mass matrix X, 



M, 



V* X 1 IT 







with X 



M 2 V2sin(3M w 
v^2 cos (3 M w /i 



where M^- is the diagonal mass matrix with the chargino masses m -± , m ~± as entries, 

"■ Xi A 2 

which are determined as the (real and positive) singular values of X. The singular value 
decomposition of X also yields results for U and V. Using the transformation matrices U and 
V, the interaction Higgsino and wino spinors iff, and W , which are two component 
Weyl spinors, can be transformed into the mass eigenstates 



%) ™th ^ 



u, 



I.) 



and ipi 



w~ 

m 



(9) 



where the ith mass eigenstate can be expressed in terms of either the Weyl spinors ipf and 



or the Dirac spinor x% ■ 

In the neutralino the neutralino mass matrix Y is symmetric, one 4x4 matrix 

is sufficient for the diagonalization 



N*YN^ = diag(m^o, m^o, m^o, m^o) 



(10) 



with 



Y 



V 



M l 
M 2 

-M z s w cos M z c w cos 

M z s w sin (3 —M z c w sin f3 



-M z s w cos (3 
M z c w cos /3 




M z s w sin (3 \ 
-M z c w sin (3 
-fi 

/ 



Mz and M w are the masses of the Z and W boson, c w = M w /M z and s w = \J\ — c^. The 
unitary 4x4 matrix N and the physical neutralino (tree-level) masses m^o (k = 1,2,3,4) 
result from a numerical Takagi factorization [52] of Y. Starting from the original bino/wi- 
no/higgsino basis, the mass eigenstates can be determined with the help of the transformation 
matrix N, 



xl 



with ^ = JV„(S°,W°, J H?, J H5 



(12) 



where ip^ denotes the two component Weyl spinor and xi ^ ne f° ur component Majorana 
spinor of the fcth neutralino field. 
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Concerning the renormalization of this sector, we implement two prescriptions that differ 
in the treatment of the complex phases. The first prescription is based on Refs. [HUES], 
while the second one is based on Ref. [2Z1EHII1HIS1]- We will emphasize the points where the 
two schemes deviate from each other. 

The following replacements of the parameters and the fields are performed according to 
the multiplicative renormalization procedure, which is formally identical for the two set-ups: 



Mi -> Mi + 6 Mi 
M 2 -> M 2 + SM 2 



1+-5Z*~ 



u~Xi 



kl 



1,2) , 
1,2) , 
(M= 1,2,3,4) , 
(M = 1,2,3,4) . 



(13) 
(14) 
(15) 

(16) 
(17) 



(19) 



It should be noted that the parameter counterterms are complex counterterms which each 
need two renormalization conditions to be fixed, (except for 5M 2 , which in scheme II is real). 
The transformation matrices are not renormalized, so that, using the notation of replacing 
a matrix by its renormalized matrix and a counterterm matrix 



X -> X 
Y -»■ Y 



6X 
SY 



with 
5X 

8Y 



5M 2 
V2 5(M w cos(3) 



5(M z s w cos p) 5(M z c w cos f3) 
\ 5(M z s w sin 0) -5(M z c w sin /3) 



v^MvKsin/?) 


<5M 2 



-5(Af z s w cos/3) 
5 (M z c w cos /3) 


— 5/i 



5 (M z s w sin /3) \ 
-5(M z c w sin /3) 





(20) 
(21) 



(22) 



(23) 



the replacements of the matrices Mf- and M^o can be expressed as 



V*5X T U f 



(24) 
(25) 



For convenience, we decompose the self-energies into left- and right-handed vector and scalar 
coefficients via 



M(P 2 )lm+^+ [^(P 2 )l m +-- M L (P 2 )lm+^ M R (P 2 )] nm ■ ( 2 6) 
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Now the coefficients of the renormalized self-energies are given by = 1, 2; k, I = 1, 2, 3, 4) 

1 



J y 



J y 



^ 2 



it 



2 L 
1 

2 L 



6Z^ + 5Z 

x x 



j ij 



5Zf_ + <5Z|! 



J y 



3fV 



£~o(p 2 ) 



X" 



J y 



3?V 



S^o(p 2 ) 



£?o(p 2 



y 

1 

fcz 2 
1 

+ 2 



^ZflM^ + iM,-5Z^ + 5M r 

-5Z L J Mt + -Mt 5Zf_ + (JMt 

2 x x 2 x x x 



5Z^o + 5Zt 



5Zt 



5Zl 



A-/ 



ki 



kl 



5M x° 



-<5Zt Mt + - 

9 X° X 9 



M^ 5Zt + 5M] 



ki 



ki 



(27) 
(28) 
(29) 

(30) 

(31) 
(32) 
(33) 

(34) 



Instead of choosing the three complex parameters Mi, M2 and \i to be independent 
parameters we impose on-shell conditions for the two chargino masses and the mass of 
the lightest neutralino (however, slightly differently in the two schemes; see below) and 
extract the expressions for the counterterms of M 1; M 2 and /j,, accordingly. It was shown 
in Ref. [371IH1E1] that for numerically stable results, one bino-, wino-, Higgsino-like particle 
should be chosen on-shell. Further, in a recent analysis [39J it was emphasized that in the 
case of the renormalization of two chargino and one neutralino mass always the most bino- 
like neutralino has to be renormalized in order to ensure numerical stability. In Ref. [IT] 
the problem of large unphysical contributions due to a non-bino like lightest neutralino is 
also discussed. In our numerical set-up, see Sec. HJ the lightest neutralino is always rather 
bino-like. On the other hand, it would be trivial to change our prescription from the lightest 
neutralino being on-shell to any other neutralino, ensuring that a neutralino with a large bino 
component is renormalized on-shell. In Ref. [3H] it was also suggested that the numerically 
most stable result is obtained via the renormalization of one chargino and two neutralinos. 
However, in our approach, this choice would lead to IR divergences, since the chargino mass 
changes (from the tree-level mass to the one-loop pole mass) by a finite shift due to the 
renormalization procedure. Using the shifted mass for the external particles and the tree- 
level mass for internal particles results in IR divergences. On the other hand, in general, 
inserting the shifted chargino mass internally yields UV divergences. Consequently, we stick 
to our choice of imposing on-shell conditions for the two charginos and one neutralino. 

As stated before, the numerical analysis is carried out using two different renormalization 
schemes. It should be noted that the differences arise in the renormalization of the param- 
eters, but the expressions for the field renormalization constants are identical. This means 
that only the renormalization of the phases of the complex parameters differ. Therefore the 
schemes are identical in the real MSSM. We will briefly describe the two schemes in the 
following, making an attempt to highlight the differences. 
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Scheme I 145,461: The on-shell conditions in this scheme read 



ReE^- (p) Xi (P) 
j u 



11 



p 2 =m 2 ■ 

Xi 



p 2 =m 2 
x l 





. 



1,2) 



(35) 
(36) 



These conditions can be rewritten in terms of six equations defining six real parameters 
or three complex ones, 



Re 
Re 



m- ± (£Uml±) + Ef-(mL) ) + EpK ± ) + E»?(mL) 



Re 
Re 



m- ± ( ^ (mi ± ) - E* (mj± ) ) - Ef (m^ ) 



,2 



X 



m *H E * o(m *;) + S *°K?. 

mjo (ijo(m|o) - Sf„(m|„; 



11 



11 










(37) 
(38) 
(39) 
(40) 



Equations (|38|) and (HOI) are related to the axial and axial-vector component of the 
renormalized self energy and therefore the l.h.s. vanishes in the case of real couplings. 
Therefore, in the rMSSM only Eqs. ( |371) and (1391) remain. It should be noted that 
since the lightest^neutralino is stable there are no absorptive contributions from its 
self energy and Re can be dropped from Eqs. (I36ll39ll40p . We retain it here in order to 
allow for these on-shell conditions to be generalized to other neutralinos. 

For the further determination of the field renormalization constants, applicable to both 
schemes, we also impose 



lim 



p 2 — >m 2 + 



lim 

p 2 —>m 2 



?)M*-(p)], 



2 

x, 



-Xi (p) 



1,2) , 



m 



ml 



(k = 1,2,3,4) 



(41) 
(42) 



which, together with Eqs. f )38|) and f l40|) . lead to the following set of equations (for 
i = 1,2; k = 1) 



Re 



E^(m? f )+Ef_K ± )) +mj ± (i£ {ml ± ) + sf_ {ml ± ) 



Xi 



Xi 



+m xf 



EfK ± ) + Ef ( m L) 



Xi 



Rc 



^-K f )-^-K f ) 



Re 



2 l S x-K ? ) 



E^K ; )j+m^ ^o K o) + E^ K? ) 

+-x? (srK ? ) + Ef'(m| ? ; 



ii 








(43) 
(44) 

(45) 
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Re 



11 



(46) 



where we have used the short-hand S'(m 2 ) = (dJ^/dp 2 
Eq. (|46p is already fulfilled due to the Majorana nature of the neutralinos 



p2 =m 2 . It should be noted that 



Inserting Eqs. (|27p - (|34l) for the renormalized self-energies in Eqs. ( 13"T|) - (T4T)]) and 
solving for [5M^-] .. and [5M^] n results in 



Im [SM^-l 
Re [5M it o] 1 

Im [5Mjo] 



1~ 
5 Re 



iRe 



E« (mL) 



-m-± Im 
2 



a? 



11 



Sf(m|o) - Sf L (m|o)] n - m,o Im[5Z,o] 



ii 



(47) 
(48) 

(49) 
(50) 



where we have used the relations ( 1441) and ( j46l) . Using Eqs. ( |22l - ( l25l) . these conditions 
lead to ESSES] 



5M l 



5M 2 



( yj | >2 v 5%o - iV*2 6M 2 + 2N* 13 N* U Sfi 

+ 2N* xl [N* 3 5(M z s w cos p) - N* 4 5(M z s w sin p)\ 
- 2N^ 2 [N* 3 5{M z c^ cos P) - N* 4 5(M z c w sin 0)} 



■x 



2 (t/i*iC/ 2 Wi^ 2 - UhU&V&VZ 

(2U; 2 V* 2 Sm- f - 2U* 12 V* 2 6m-± 

+ (Wi " ^1^)^22 S(V2M W sin/3) 
+ Ut 2 U; 2{ V^-V^)6(V2M wC osP) 
1 

2 (Hjisjwfo - u* 12 u* 21 v? 2 v 2 \ 



■x 



2^^ 5m-± - 2U* 21 V 2 \ 5m^± 



+ (CM " cos /3) 

+ " n*^2* 2 ) 5(^2% sin /3) 



(51) 



(52) 



(53) 



where, combining Eqs. ( 14~7l) - ( 1481) and Eqs. (1491) - (l50l) . we introduced the short-hand 
notation 



dm^o 

Al 



-Re 



Xr 



11 



(54) 
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Sm xt = l Re 



(55) 



Here we have used already Eqs. ( 173!) and (176"]) . see below, which fix the expressions 
for the imaginary parts of the chargino and neutralino field renormalization constants 
5Z^_ and <5Z^o, respectively. 

Scheme II [371138] : Here on the other hand, it is the real part of the corrections to 
the self-energies of the on-shell particles that are required to vanish, 



( Re ±x- (p) Xi (p) 



ReS^o(p) xiip) 



ii 



p 2 =m 2 , 



9 9 



o 

. 



1,2) 



(56) 
(57) 



It should be noted that for the derivation of the field renormalization constants, in 
order to ensure that the on-shell propagator has only a scalar and vector part, we 
impose the additional conditions 



U L (m 2 



Xj ■ 



33 



3.1 



U R (m 2 



(58) 
(59) 



33 



with Xj denoting either a chargino (with j = 1,2) or a neutralino (with j = 1, . . . , 4). 
The first equation, relating the vector coefficients, is automatically satisfied in the 
cMSSM for both charginos and neutralinos. This means that the conditions (jUJ) 
and (146!) applied in scheme I, describing the axial and axial-vector components of the 
renormalized self-energy, are also satisfied in scheme II. It should be noted that here, 
however, we have dropped the R0. On expanding Eqs. (15"6"|) and (foTl) . in analogy to 
Eqs. ([37]) and we find 



Re 



Re 



+ Sf(^) + Sf(^) 

+ ^k ; ) + ^k ? )" 



0, 
. 



11 



(60) 
(61) 



We are left with three on-shell conditions for scheme II, Eqs. ( l60|) and ( jBTT) . Therefore 
one can only fix the renormalization constant^! 5|Mi|, 5M2 (note that Mi is chosen to 
be real) and <5|/z|, but not 5</?Mi or where 5Mi = (5|Mi|e l</3M i and Sfi = 5\ix\e 1 ^. 
This, however, is not a problem, as it turns out that the phases of M\ and /i are UV 
finite at one- loop (see e.g. Ref. [37]), and need not be renormalized at all, i.e. they 
can be set to zero, = 5(f>^ = 0. It should be noted that 5\/jl\, S\Mi\, SM 2 are 

^ee the discussion on absorptive contributions at the end of this subsection. 

2 Here we adopt the notation S for scheme II to distinguish the renormalization constants from those in 
scheme I. 
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related to 5M^- and 5M^o via Eqs. ( 122|) - ( 125]) . Expressions for these renormalization 
constants can then easily be obtained by inserting Eqs. (J27J) - (13"4"|) for the renormalized 
self-energies in Eqs. ( 16"Uj) - flBTj) . resulting in (for details see Ref. [37JI3B1III1I2I]) 



Re( e -^M 1 N 2^ s 



x 



5M 2 



[Re (U n V n ) Re (e-^U 22 V 22 ) - Re (e-^[/ 12 V 12 ) Re (U 21 V 2l )] N x 
+ [2 Re (e'^NuNu) Re (*7nVn) + Re (JV? 2 ) Re (e^?7 12 Fi 2 )] C 2 
- [Re (A^) Re ( e -^[/ 22 V 22 ) + 2 Re (e^N w N u ) Re (17 21 V 21 )] d 

i [Re (e-^?7 12 y 12 )C 2 - Re (e- i ^U 22 V 22 )C 1 ] , 

--^[Re^n^iOd-Re^^Od] , 



(62) 
(63) 
(64) 



where we use the abbreviations 



Q = Re 



25(M w sin/3) Re (U i2 V a ) - 25(M W sin/3) Re (U a V i2 ), 



(65) 



iV, = Re 



m,o[E^ (m|o) + SfoKo)] + Sf L (m|o) + Sf R K ) 



y v u v""v ( v 1 "v u v""v u yj 1 "y u v""v u y 1 "v u v — v w 

■^-Z A A-i A Aj A Aj A Aj 

+ 45(M z s w cos (3) Re (iV a ^ 3 ) - 45(M z c w cos (3) Re (iV i2 iV i3 ) 
- 45(M z s w sin /3) Re (NuNa) + 45(M zCw sin /?) Re (N i2 N u ) , 

5 = 2[Re {U 21 V 21 ) Re (e^l/iaVy - Re (U n V n ) Re (e^E/aaV^)] 



(66) 
(67) 



5M^- and 5M^o are simply obtained by the replacements 5 Mi, SM 2 , and 5fi by, 
respectively SM U SM 2 , and 5// in Eqs. fl20|) - ( 125|) . 

The following discussion of the neutralino mass shifts and the field renormalization is appli- 
cable to both schemes. 

Since the chareino masses m~±,m~± and the lightest neutralino mass m^o have been 

° Xl ' X2 Xl 

chosen as independent parameters, the one-loop masses of the heavier neutralinos Xi = 
2,3,4) are obtained from the tree-level ones via the shifts 



A.J Aj 



Am.o = Re { m*o ( S^(mL) + S"„(mL) ) + Eff (mi,) + £^(m~o 



Xi • 



A^ A^ 



Where necessary we distinguish the tree-level mass rriyO from the on-shell mass 

Aj 



m i o = m^o + Am^o 

Ai Ai Xi 



(6£ 



(69) 



We use m^o for all externally appearing neutralino masses, which includes the (on-shell) 
momentum in the employed neutralino self-energies. In order to yield UV-finite results we use 
the tree-level values m^o for all internally appearing neutralino masses in loop calculations. 
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Eqs. (143!) and (j4"5l) define the real part of the diagonal field renormalization constants of 
the chargino fields and of the lightest neutralino field. By extending Eq. (l4"o*j) to apply also 
for k = 2,3,4, we can generalize the result for the diagonal field renormalization constants 
of the lightest neutralino to the other neutralino fields. 

The imaginary parts of the diagonal field renormalization constants are still undefined. 
However, these can be obtained using Eqs. f|4*8|) and (1501) . where the latter is generalized 
to include k = 2,3,4. Now in scheme I (or scheme II), for the charginos and the lightest 

neutralino, Eqs. (jUJ) and fl50l define the imaginary parts of [5M^±] .. (or 5M^± ) {i = 

1,2), and [M ? ] n (or 5M ^0 ) in terms of the imaginary part of the field renormalization 

constants. Therefore these are simply set to zero (see below Eqs. (1731 and fl7"6l) ). which is 
possible as all divergences are absorbed by other counterterms. 

The off-diagonal field renormalization constants are fixed by the condition that 



ReE^- (p) . . Xi ip) 



ReE^o(p) Xi(p) 

J kl 



p 2 =m 2 , 



9 9 

p z =m n 



= 1,2), 

(fc,/ = 1,2,3,4) 



(70) 
(71) 



Finally, this yields for the field renormalization constants [53J (where we now make the 
correct dependence on tree-level and on- shell masses explicit), 



Re 



Im 



5Z x 



L/R 



5Z 



L/R 



5Z 



L/R 



Re 



Im 



-Re 



E L/R (m 2 ± ) 



V 



:mJ±)+^(my)+,/< ;= (!: 



SL' 











t 1 




- ii 







m 



- E 



Xi y X 

Im SM 



X 



Si 



2 



2 



Re 



m 2 ~±Yi^_ R (m 2 ~±) + m-±m^±H~_ L (m~±) 

v. v \ v. / y.- y_- y V v. / 



. ^SL/SR, 2 \ , ^SR/SL, 2 

+ m-±E-_ (m-±) + m~±E~_ (m_ r 



m-± 5M 

Xj/i 





kk 


-Re 














1 




kk 


Xk 



4°/ 



Xk ' 



+ m^ ( E^o(mto) + E$)(mSo) ) + m^o ( E^ (mi ) + E5f (m^ 



Xfc ' X 
k=l, Si 

J kk 



(72) 

j 

ii 

(73) 
(74) 



(75) 



A' A- 



(76) 



5Z- 



J fcz 



m-n — m 2 - 

fc X; 



Re 



«4? E y°(^y9) + ^™# S ?o(m?p) 



+ m^oS^(mto) + rn^o E£o (mto) - m^SM^ - m^o5Ml 



(77) 



within scheme I (Si). Making the replacements 5M^o — y 5M^o and SNL^- — >■ 5M^- for 
scheme II, Eqs. fl7"3l) and f!76|) no longer vanish. 
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Contributions to the partial decay widths can arise from the product of the imaginary 
parts of the loop-functions (absorptive contributions) of the self-energy type contributions in 
the external legs and the imaginary parts of complex couplings entering the decay vertex or 
the self-energies. It is possible to combine these additional contributions with the field renor- 
malization constants in a single "Z factor" , Z, see e.g. Ref. |l5lll6] and references therein. In 
our notation they read (unbarred for an incoming neutralino or a negative chargino, barred 
for an outgoing neutralino or negative chargino, and not making the difference between 
scheme I and II explicit), 



8Z. 



L/R 



E L/R (m 2 . ± : 

X Xi 



+ m xt 



E^'fmD + Ef? K ± ) 



Xi 



(7f 



8Z. 



L/R 



L/R 



kk 



x° 



L/R 



kl 



± 



2m xt 



+ 5Mt_ 



9 

m + 

Xi 



2 

m + 

Xi 



2 ^L R, 2 \ i L t 2 \ 

m~±h~L \m.± ) + m~±m~±h-_ (m~±) 
xf X v xf J Xi Xj X v xf> 



+ m-±zZ S ; L J SR (m 2 ± ) + m-±Zf R/SL (m 2 ± ) - m-± 8M 9 - - m~± 5M\_ 



Xi X 



Xi/j 



(79) 



30) 



+ m 2 - 

Xk 



± 



Eg(rh|o) + Sg(m|o; 



+ m xl 



(s-'(^o) + S-'(^o)) 



kk 



2m^o . 

Xk 



m 2 

Xk 



m 2 
x° 



m z ~oX L/ R (mi ) + m f om f oS^(m-o) 

X; X° V Xi> X k Xi x K Xi> 



R/L . 



6M h 



J kk 



(81) 



+ m^oS^ L/5i? (m|o)+m^oS 



SR/SL /~2 s 
o (m-o, 



m^o <5MvO — m^o 5Mt Q 

Xfc/i x x l/k x u 



L/R 



5Z 



L/R 



E L J R (ml ± ] 

X Xi 



(82) 



+ m xt 



2 's£fm 2 _ ± 



2m xt 



1 r ^K i± )-Sf«(m 2 i± )-^- + 5Mt 



2 2 

m:± — mt± 

Xj X 



m lt ^ R ^ m \t } + m xt m xf ^ ± ■ 



Xi 



f)7 L/R 
oz, x° 


kk 


b? R/L 
°*X° 




8Z L ' R 
x" 


kl 


8Z R i L 

x u 





(83) 
i 

(84) 
(85) 
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within scheme I, and with iM^o — > <5M^o and <5M^,- — y SlSA^,- for scheme II. The char- 

gino/neutralino Z factors obey KeSZ^ R = [Re8Z^' R ]^ = [5Z^ R ]\ which is exactly the 
case without absorptive contributions. The Eqs. fISll) and fISol) hold due to the Majorana 
character of the neutralinos. We will use these Z factors rather than the field renormalization 
constants defined in Eqs. ( !72|) to (1771) in the following numerical analysis. 

Special care has to be taken in the regions of the cMSSM parameter space where the 
gaugino-Higgsino mixing in the chargino sector is maximal, i.e. where ~ M 2 . Here 5M 2 
(see Eq. §2$) and 5fi (see Eq. (E3}) diverge as (U^U^V^V^ - U* 2 U; x _1 for scheme I 
or 2[Re (U 21 V 21 ) Re {e'^U^V^) - Re {U-nVu) Re (e^U^Va)] for scheme II, and the loop 
calculation does not yield a reliable resulto. It should be noted that the singularity arises in 
both schemes for = M 2 . These kind of divergences were also discussed in Ref. [3~y [l39lHI] . 

Our two renormalization schemes differ in the treatment of complex contributions in the 
chargino and neutralino sector. The first difference is that in scheme I, in the derivation of 
the renormalization conditions 5M 2 is allowed to be complex. Therefore six real conditions 
must be imposed in order to renormalize the mass matrices X and Y, Eqs. ([H]) and ffTT]) . as 
opposed to the five in scheme II. In addition, in scheme I, the renormalization of the phases 
is obtained by imposing that the imaginary parts of the relevant diagonal renormalization 
constants Eqs. (!73j) and (1761) vanish. In scheme II the phases are not renormalized as they 
are found to be UV finite. 

Scheme I is based on the idea that the absorptive contributions should not enter the 
renormalization procedure. Therefore it requires that the absorptive contributions are not 
included in the on-shell mass renormalization conditions Eqs. fl37|) and fl39l) . While the 
phases of the complex counterterms of the mass matrix parameters are also found to be 
UV finite, no condition is imposed on them. Notice that it is always possible to rephase 
the parameters of the chagino/neutralino sector since only the relative phases are physically 
relevant. Scheme II is supported by the argument that the phases are not renormalized, 
i.e. they are the same at tree level and loop level. Therefore it is clear which of the funda- 
mental parameters are chosen to be real and which complex, i.e. in scheme II M 2 is fixed 
to be real while in scheme I its counterterm is allowed to be complex. As stated earlier, the 
results in both schemes should agree in the case of real parameters. In the complex case, we 
expect the differences to be very small, as they are of higher order. 

2.2 The lepton/slepton sector of the cMSSM 

For the discussion of the one-loop contributions to the decay channels in Eqs. ([5]) and ([6]) 
a description of the scalar lepton (£) and neutrino {pi) sector as well as their fermionic SM 
partners is needed (we assume no generation mixing and discuss the case for one generation 
only). The bilinear part of the I and vi Lagrangian, 

= - (4 4) (|) - (#) M u - e (u t ) , (86) 

3 Similar divergences appearing in the on-shell renormalization in the sbottom sector, occurring for "max- 
imal sbottom mixing", have been observed and discussed in Refs. [S5 1 H5 1 I5B]. 
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contains the slepton and sneutrino mass matrices and M^, given by 



M 
M 



Mf + mj + M 2 z c 2/3 (I$ -Q t sl 



m t Xi 



"I 



Mj L +llc w M 2 z 



Mf + mf + M z c 2 pQ e sl 



57) 



with 



Xi = An — fi* tan f3 



(89) 



Mji L and Mj, are the soft SUSY-breaking mass parameters, where M| is equal for all 



members of an SU(2)l doublet, me and Qi are, respectively, the mass and the charge of the 
corresponding lepton, If, denotes the isospin of £/v, and Ae is the trilinear soft-breaking 
parameter. We use the short-hand notations c x = cos(a;), s x = sin (a;). The mass matrix 
can be diagonalized with the help of a unitary transformation U^, 



D,-=U,-M,-U 




Ui 



Us u, 



(.22 



(90) 



The mass eigenvalues depend only on \Xi 
ordered, i.e. < : 



The scalar lepton masses will always be mass 



h,2 2 



\ (M\ L + MfJ + m 2 + \tfc w Ml 



=F 



Mf - Mf + M z c 2p {If - 2Q l8 %) + Am 2 \X e \ 2 , 



m% = Mf + llc w M z . 



(91) 



(92) 



A detailed description of the renormalization of this sector can be found in Refs. jlHlHHJET] . 
Here we just briefly review the restoration of the SU(2)l relation for the renormalized slep- 
ton mass parameters. Since our on-shell approach results in an independent renormalization 
of the charged sleptons and of the scalar neutrino, we need to restore the SU{2) i relation 
at one-loop to avoid problems concerning UV- and IR-finiteness as discussed in detail in 
Ref . [56J . This is achieved via a shift in the M| parameter entering the I mass matrix (see 
also Refs. [l^|l6j[57J[58] ) • Requiring the SU{2) L relation to be valid at the loop level induces 
the following shift in Mf (£) 



with 



Mf T (I) = Mf T [y t ) + 6 Mf r {y t ) - bMffi) 



(93) 



SMffi) 



SMf L (u £ ) 



+ M 2 z c 2 p QtSsl 



2 r) *~ 2 (II - Q(sl){c w 6M 2 Z + M\ 5c w ) , 



6m%-I*{c 2 p6M z + M z 5c 2 p) 



(94) 
(95) 



14 



Such shifts however mean that both slepton masses are no longer on-shell. An additional 
shift in Mj restores at least one slepton mass to be on-shell. 



M W = TT^VT Tl'^L ^Tf ? n 2 , 2 ~ m * ~ M * ^ Q< s ^ + < ■ ( 96 ) 

Mf L (£) + mj + Mf c 2/3 (Jf - Q t s%) - m? 



A "natural" choice is to preserve the character of the sleptons in the renormalization process, 
and this additional shift relates the mass of the chosen slepton to the slepton parameter M^ r . 
As mi x < mi % (see above) and < M~ for both scenarios considered later (see Tab. [I]), we 
choose to insert m| into Eq. (}96l) and recover its original value from the re-diagonalization 
after applying this shift. 

The renormalization of the quark/squark sector is described in detail in Ref. [16]. As for 
the slepton sector, the restoration of the SU(2)l relation in the squark sector leads to shifts 
analogous to Eqs. (|9~5j) and (|9~6"j) in the left and right squark mass parameters, respectively, 
Mg L and Mg R . In our subsequent numerical analysis these shifts have been included for the 
results calculated in scheme I, but their effect is found to be negligible for our choice of 
parameters. 

2.3 The Higgs boson sector of the cMSSM 

The two Higgs doublets of the cMSSM are decomposed in the following way, 



H 



12 



'1 



= (u 21 ) =e*[ i Z m ■ J • ( 97 ) 



22 



Besides the vacuum expectation values Vi and v 2 , in Eq. (j97|) a possible new phase £ between 
the two Higgs doublets is introduced. After a rotation to the physical fields one obtains for 
the terms linear and bilinear in the fields, 



V H = . . . + T h h + T H H + T A A 

fh\ 
H 

A 

\GJ 



- ^ (h, H, A, G) Mf^ AG 



+ {H\G+)Mf£ G± ^+... , (9? 



where the tree-level masses are denoted as nth, ran, rriA, tug, M h ±, m G ±. With the help 
of a Peccei-Quinn transformation [59] [i and the complex soft SUSY-breaking parameters in 
the Higgs sector can be redefined [50] such that the complex phases vanish at tree-level. 
Including higher-order corrections the three neutral Higgs bosons can mix [5| [TUllTi] l5T)] . 

(h, H, A) — > (h u h 2 ,h 3 ) , (99) 

where we define the loop corrected masses according to 

M hl < M h2 < M h3 . (100) 
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Details about the renormalization and the Z factors, ensuring the on-shell properties of 
external Higgs bosons can be found in Refs. [4"6ll3U] . For the renormalization of tan /3 and the 
Higgs field renormalization the DR scheme is chosen [321130]. This leads to the introduction 
of the scale which will be fixed later to the mass of the decaying particle. 

3 Calculation of loop diagrams 

In this section we give some details about the calculation of the higher-order corrections 
to the neutralino decays. Sample diagrams are shown in Figs. [TJ [61 Here the generic 
internal particles are labeled as follows: F denotes a SM fermion, chargino, or neutralino, 
S denotes a sfermion or a Higgs and V denotes a 7, Z, or W ± . Concerning the diagrams 
for decays into charged particles, although we only show diagrams for decays into Xi or 
leptons, we also include the corresponding diagrams for the decays into xt or antileptons. 
Not shown are the diagrams for real (hard or soft) photon radiation. These are obtained 
from the corresponding tree-level diagrams by attaching a photon to the electrically charged 
external particles. It should be noted that the expressions for the tree-level diagrams are 
given explicitly in Appendix |A] Counterterm diagrams are also not shown, but and these 
can be obtained from the corresponding tree-level diagrams, replacing the tree-level vertex 
by the counterterm vertex. 

Internally appearing Higgs bosons do not receive higher-order corrections in their masses 
or couplings, which would correspond to effects beyond one-loop. Furthermore, we found that 
using loop corrected Higgs boson masses and couplings for the internal Higgs bosons leads 
to a divergent result. For external Higgs bosons, as mentioned in Sec. 12.31 the appropriate 
Z factors are applied, following the prescription of Ref . [30] . For the numerical analysis, these 
factors, as well as the loop corrected masses for external Higgs bosons, are obtained from 
FeynHiggs-2.9.0. Diagrams with a gauge boson/ Goldstone-Higgs self-energy contribution 
on the external Higgs boson leg, absent from the Higgs Z factors [30], are also required for 
the decays x° — > x?^fc> (i — 2,3,4; j < i; k — 1,2,3), Fig. (H with a Z/G-h k transition, 
and similarly for the decay x! -> Xj H+ (» = 2 > 3 > 4 5 3 = M), Fig. El with a W + /G + -H + 
transition. On the other hand, Goldstone-Higgs/gauge boson self-energy corrections for the 
neutralino decay to a chargino/neutralino and a gauge boson xl ~ * X^Z (i = 2,3,4; j < i), 
or x°i XjW + , (i = 2,3,4; j = 1,2), can be neglected, as these vanish on mass shell, i.e. 
for p 2 = M| (p 2 = M$y) due to e -p — 0, where p denotes the external momentum and e the 
polarization vector of the gauge boson. 

The diagrams and corresponding amplitudes have been obtained with FeynArts [61]. The 
model file for calculations in scheme I, including the MSSM counterterms, is described in 
more detail in Ref. |46j , and the model file used for calculations in scheme II is based on that 
discussed in Refs. [2Z1I33]- The further evaluation has been performed with FormCalc (and 
LoopTools) [62J. As regularization scheme for the UV-divergences we have used constrained 
differential renormalization [63J , which has been shown to be equivalent to dimensional reduc- 
tion [64J at the one-loop level [62]. Thus the employed regularization preserves SUSY [631166]. 
All UV-divergences cancel in the final result. 

The IR-divergences from diagrams with an internal photon have to cancel with the ones 
from the corresponding real soft radiation, where we have included the soft photon con- 
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Figure 1: Generic Feynman diagrams for the decay — > x°jhk {i — 2,3,4; j < i; k = 
1,2,3). F can be a SM fermion, chargino, or neutralino; 5* can be a sfermion or a Higgs 
boson; V can be a 7, Z, or W ± . Not shown are the diagrams with a Z—hk or G-h^ transition 
contribution on the external Higgs boson leg. 




Z Z Z 




Z 



Figure 2: Generic Feynman diagrams for the decay %° — > \^-Z (i = 2,3,4; j < i). F can 
be a SM fermion, chargino, or neutralino; S can be a sfermion or a Higgs boson; V can be 
a 7, Z, or W ± . 

tribution following the description given in Ref. [67]. The IR- divergences arising from the 
diagrams involving a 7 are regularized by introducing a finite photon mass, A. All IR- 
divergences, i.e. all divergences in the limit A — > 0, cancel to all orders once virtual and real 
diagrams for one decay channel are added. The only exception are the decays X2 3 4 ~ ^ xf 
The shift to the neutralino on-shell masses via Eq. (1681) results in an IR divergence at the 
two-loop level, i.e. here we find a cancellation of the divergences "only" at the one-loop level, 
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H+ H+ H+ H+ 




Figure 3: Generic Feynman diagrams for the decay Xi ~^ XjH + (i = 2,3,4; j = 1,2). 
F can be a SM fermion, chargino, or neutralino; S can be a sfermion or a Higgs boson; V 
can be a 7, Z, or W . Not shown are the diagrams with a W + -H + or G + -H + transition 
contribution on the external Higgs boson leg. 




WWW 




W W 



Figure 4: Generic Feynman diagrams for the decay Xi ~^ XjW + (i = 2,3,4; j = 1,2). F 
can be a SM fermion, chargino, or neutralino; S can be a sfermion or a Higgs boson; V can 
be a 7, Z, or W ± . 

as required for our one-loop calculation. The remaining two-loop IR divergence could be 
eliminated by a symmetry restoring counterterm in the X234 — xfW~ vertex, similar to 
the evaluation of the decay t 2 — > bi^W + in Ref. |46]. We have furthermore checked that our 
result does not depend on AE defining the energy cut that separates the soft from the hard 
radiation. Our numerical results have been obtained for AE = 10 -5 x m^o for all channels^ 

4 The larger cut is necessary to obtain a better convergence of the integration over the three body phase 
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i- t~ r 




i- t~ 



Figure 5: Generic Feynman diagrams for the decay x° — > i~ l\ {% = 2, 3, 4; i = e, fx,r; k = 
1, 2). F can be a SM fermion, chargino, or neutralino; S can be a sfermion or a Higgs boson; 
V can be a 7, Z, or W ± . 



v\ v\ v\ v\ 




Figure 6: Generic Feynman diagrams for the decay x° — > vt.v\ (i = 2,3,4; i = e, fx, r). F 
can be a SM fermion, chargino, or neutralino; S can be a sfermion or a Higgs boson; V can 
be a 7, Z, or W ± . 

4 Numerical analysis 

In this section we will first introduce and motivate the scenarios studied, discussing the 
current experimental constraints considered, then introduce the observables calculated, and 
finally present our results for each of the decay channels of the heavier neutralino (%° — > xy) 
as a function of y?Mi • 

space. The contribution from nearly collinear photons (along the direction of the electron) leads to numerical 
instabilities in the integration. 
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M H ± 












^& 




4r 


20 


160 


600 


350 


300 


310 


400 


1300 


1100 


2000 



Table 1: MSSM parameters for the initial numerical investigation; all mass parameters are 
in GeV. Mi, M 2 and /x are chosen such that the values for m~± and m~± and Eq. (11011) are 
fulfilled (see text). 

As stated earlier, we present our results in two scenarios. In both, the absolute value of Mi 
(see above) is fixed via the GUT relation (with |M 2 | = M 2 ) 

| Mi | = ^tan 2 # w M 2 « \m 2 . (101) 

For the numerical analysis we obtain M 2 and /i from the fixed chargino masses m y.f 2 , and 
| M x | via Eq. (11011) . leaving ip^ as a free parameter. Our two scenarios arise due to the 
ambiguity in calculating \x and M 2 from Tn^± • This ambiguity can be resolved by choosing 
an addition condition, /j > M 2 or // < M 2 . The first choice, denoted by <S/i, results in a 
Higgsino-like x° while the second choice, denoted by S g , results in a gaugino-like xl- 

The values of the parameters for these scenarios are given in Tab. (U where, in analogy to 
the slepton parameters M^ , Mjj R and A t defined in Sec. 12.21 for the sleptons, Mg L and Mg R 
are the left- and right-handed soft SUSY-breaking mass parameters and A q is the trilinear 
soft-breaking parameter for the squarks. These are chosen such that most decay modes are 
open simultaneously to permit an analysis of as many channels as possible. Only decays into 
the heavier chargino xfi 1X1 general degenerate with the heavier neutralino, and the decay 
channels Xz^k/Z, (k = 1, 2, 3), in Sh are kinematically closed. We also ensure that the 

scenarios are consistent with the MSSM Higgs boson searches at LEP [3], Tevatron [68] and 
LHC [69] • The light Higgs mass scale together with the value of tan/3 = 20 are in potential 
conflict with the recent MSSM Higgs search results, which, however, have only been obtained 
in the m™ ax scenario [70] . We stick to our parameter combination to facilitate the numerical 
analysis with all decay channels involving Higgs bosons being open simultaneously. On 
the other hand the recent discovery of the (lightest) Higgs boson at the LHC [I] allows for 
tan (3 > 9 as given in Tab. [T]@ Furthermore, the following exclusion limits for neutralinos [71] 
hold in our numerical scenarios: 

m^o > 46 GeV, m^o > 62 GeV, m^o > 100 GeV, m^o > 116 GeV . (102) 

It should be noted that the limit for m^o arises solely when Eq. (11 011) is assumed to hold. 
In the absence of this condition, no limit on a light neutralino mass exists, see Ref. [72] and 
references therein. 

The most restrictive experimental constraints on the phase of (pM 1 arise due to the 
bounds on the electric dipole moments (EDM's) of the neutron d n , mercury rf Hg , and Thal- 
lium d Ti [75H75] FI Using CPsuperH2.2 [SD], we have calculated these EDM's, and find that 

5 The Higgs mass in the allowed range can be obtained varying the squark trilinear coupling A q . 
6 In addition, the heavy quarks [75], the electron [771[7S] and the deuteron [75] EDM's should be taken 
into account. 
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Scenario 




m xt 








A-l 


/i 


M 2 


Mi 




600.0 


350.0 


600.0 


364.2 


359.6 


267.2 


362.1 


581.8 


277.7 


s h 


600.0 


350.0 


600.1 


586.2 


349.9 


171.4 


581.8 


362.1 


172.8 



Table 2: The chargino and neutralino masses in the scenarios S g and Sh- We also show 
the values for the "derived" parameters Mi, M 2 and /x. All mass parameters are in GeV, 
rounded to 0.1 GeV to show the size of small mass differences, which can determine whether 
a certain decay channel is kinematically closed or open. 



in the scenarios studied the bounds due to the Thallium EDM are the most constraining, 
ipMi ^ vr/100. This is mainly because, in order to keep all channels open the selectrons 
in our scenarios are light. On increasing the masses of the lower generation of sleptons to 
1.2 TeV (f Mx is unconstrained. It should be noted that such a change would only affect our 
results such that the decays to the lower generation slepton channels are no longer open: the 
loop corrections to the other decays are largely independent of these masses. 

The chargino and neutralino masses for <£>Mi — are shown in Tab. [2], while the Higgs and 
slepton masses are shown in Tab. [3j Here h 2 corresponds to the pure CP-odd Higgs boson. 
For (fiMi — 90° (i.e. the maximal CP-violation possible in our numerical analysis) we find the 
same Higgs boson masses within the precision of Tab. El In this case /i 2 receives a very small 
CP-even admixture of < 0.003% in both scenarios, while h\ and h 3 remain correspondingly 
a nearly pure CP-even state. The masses m^o are chosen such that the neutralinos would be 
copiously produced in SUSY cascades at the LHC. Furthermore, the production of XiXp f° r 
i = 1,2 and j = 2,3,4, at the ILC(1000), i.e. with yfs = 1000 GeV, via e + e' -> jfixj will 
be possible, where unpolarized tree-level cross sections in the scenarios S g and Sh are shown 
in Tab. HI All the subsequent decay modes (EI) - ([6]) would be (in principle) open, and the 
clean environment would permit a detailed study of neutralino decays [15JEI6] . Higher-order 
corrections to the production cross sections would change the values in Tab. 0] by up to 
0(10%) [81] . and choosing appropriate polarized beams, could enhance the cross-sections by 
a factor 2 to 3. The accuracy of the relative branching ratio Eq. (11051) at the ILC would be 
close to the statistical uncertainty, and from the high-luminosity running of the ILC(1000), 
a determination of the branching ratios at the percent level might be achievable. We have 



Scenario 




mfi 2 


m fl 




m ^ 


mt T 


M H ± 


m hl 


m h2 


m h3 




303.2 


313.1 


287.3 


328.0 


293.0 


293.0 


160.0 


125.8 


137.2 


140.3 


s h 


302.9 


313.3 


273.7 


339.5 


293.0 


293.0 


160.0 


125.8 


137.4 


140.3 



Table 3: The slepton and Higgs masses in the scenarios S g and Sh- The selectron and electron 
sneutrino masses are equal to those of the corresponding smuon and muon sneutrino up to 
a few tenths of GeV. All masses are in GeV, rounded to 0.1 GeV. 
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Seen. 


process 


o-o.q [fb] 


fpoiffb] 


stat. prec.0,0 


stat. precp i 




e+e" -> xlXi 


5.2 


15.0 


4% 


3% 


s h 


e+e- xlXi 


1.0 


1.6 


10% 


8% 


s 9 


e+e- -> X4X2 


1.0 


2.9 


10% 


6% 


s h 


e+e" -> X4X2 


0.4 


1.1 


16% 


10% 


s 9 


e+e" -> X4X3 


0.5 


0.4 


14% 


16% 


s 9 


e + e" ->■ E 


6.7 


18.3 


4% 


2% 


s h 


e + e" ->■ E 


1.4 


2.7 


8% 


6% 


s 9 


e + e" E 


39.9 


71.8 


2% 


1% 


s h 


e + e" E 


3.2 


6.1 


6% 


4% 


s 9 


e + e" E jgxS 


13.5 


11.0 


3% 


3% 


s h 


e + e" E 


46.2 


132.4 


1% 


1% 



Table 4: Neutralino production cross sections at the ILC 1000. Here o^o denotes the cross 
section for unpolarized beams, while <7 po i denotes that with electron and positron polarization 
—80% and +60%, respectively. The sums of XiXp f° r * = 2,3,4, are performed over j < i. 
The two right-most columns show the statistical precision for a (hypothetical) branching 
ratio of 10% assuming an integrated luminosity of 1 ab -1 , rounded to 1%. 



calculated the decay width at tree- level ("tree") and at the one- loop level ("full"), including 
all one-loop contributions as described in Sec. El and in addition the relative size of this 
one-loop correction via, 

"pfull "ptree 

r tree = r tree ( ~0 ^ xy) ? pfull = p fu ll ( ~0 ^ xy) ? Ar/r = ^ . (103) 

In the figures below we show the absolute value of the various decay widths, T(xl ~ > x j) on 
the left and the relative correction from the full one- loop contributions on the right. The 
total decay width is defined as the sum of all kinematically open two-body decay widths, 

rf a f = J2 r trce (x° -> xy) , = £ r full (x° -> xy ) . (104) 

xy xy 

The absolute and relative changes of the branching ratios are defined as follows, 

DD u,o _ r—tx? -> xy) M _ r M '(x? -> xy) ABR _ BR MI - BR"" 

BR = fff ' = TS ' "BR" = BR'- ' (1 ° 5) 
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The last quantity is crucial in order to analyze the impact of the one-loop corrections on 
the phenomenology at the LHC and the ILC, see below. Since decays to a light Higgs 
and the LSP are of particular importance, for the decay x\ ~~ >* Xi^i, we also show, in the 
lower panels, the branching ratio BR(%4 — > xy) (left) and the relative size of the one-loop 
correction (right). The corresponding branching ratios for the other channels can be inferred 
from these plots. 

In order to distinguish the results evaluated in the two schemes we denote those of 
scheme II with a 'tilde', i.e. S g and Sh for the scenarios S g and Sh, respectively. It should 
be noted that the tree-level results obtained for the two schemes fully agree, as our two 
renormalization schemes differ only in the treatment of complex parts. The difference of 
one-loop results for real parameters is negligible, and is only due to a different handling of 
the corrections in the squark sector, which is not highly relevant for the electroweak decays. 
Therefore the results for both schemes are only shown on the right panels for the relative 
corrections. In Sec. I4.6l we will summarize and discuss the differences between these schemes, 
highlighting those channels where the deviations are largest. 

The numerical results we show in this section are of course dependent on the choice of 
the MSSM parameters. Nevertheless, they give an idea of the relevance of the full one-loop 
corrections. Decay widths (and their respective one-loop corrections) that may appear to be 
unobservable due to the small size of their BR, could become important if other channels are 
kinematically forbidden. Consequently, the one-loop corrections to all channels are evaluated 
analytically, but in the numerical analysis we only show the channels that are kinematically 
open in our numerical scenarios, except for the decays into leptons of the first two families, 
which are closely related to the decays into third family leptons. 

4.1 Decays into charged Higgs and W bosons 

We start our numerical analysis with the decays x\ ~~ X\H + , presented in Fig. [7J The 
partial decay width for the charge conjugated process xl ~~ >* X\H~ can be obtained by 
taking the charge conjugate of all the couplings. Since in our analysis only M\ is complex 
this is obtained with the transformation </?Mj 2tt — ifM x - This argument will also be valid 
for all the decays into VT-bosons and lepton-slepton pairs described below. All these decays 
have also been computed in order to evaluate the total decay width at one-loop level. 

This channel yields decay widths of around 0.7 GeV in both scenarios, corresponding to 
BR's of ~ 4.5% for S g and ~ 11% for Sh- The tree-level partial widths are almost equal in 
both scenarios due to the symmetry of the Higgs-gaugino-Higgsino couplings under <->■ \i 
exchange, as well as to the similar phase space. The relative corrections, shown in the right 
plot, are of the order of a few percent. The mild dependence on (p^ for this process is 
due to our choice of parameters, in particular the GUT relation on the gaugino parameters 
Mi and M2, which leads to a weak dependence of the heavier neutralinos on M±. The dips 
best visible in the right panel are due to the Xi XiW~ threshold, at <£mi — 125° and 
ip Ml ~ 235°. Due to CPT-invariance the masses are invariant under ip Ml — » — resulting 
in mirrored threshold effects. These effects are further discussed in Sec. 14.21 for the decays 
into the second lightest neutralino. 

The results for the two schemes are shown for the relative corrections in the right panel. 
The agreement between the relative corrections is at the level of 10 -5 for both scenarios, 
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Figure 7: r(x° — >■ X\H + ). Tree- level ("tree") and full one- loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [H with varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 



see Tab. El and therefore cannot be visibly distinguished here. It should be noted that the 
differences between the schemes are particularly small due to mild dependence of the tree 
level decay width on <fMi- As the schemes are identical in the real case, only decays with a 
stronger dependence on the phase ifMt are sensitive to the differences between them. 

Next we analyze the decays x\ — ^ XiW + shown in Figs. [HJ The general behavior of the 
decays into W + is very similar to those into H + discussed above. This decay yields decay 
widths around ~ 1 GeV in both scenarios, corresponding to BR's of ~ 7% and ~ 17% for, 
respectively, S g and Sh, with a mild dependence on ifM x - The one-loop effects are found to 
be, respectively, of ~ —8% and ~ —4% and the same thresholds as in the previous process 
can be observed as dips in the right panels. The agreement between the relative corrections 
is at the level of 1CT 5 for S g and 10 -4 for Sh, see Tab. EJ and therefore is again too high to 
be observed. 

4.2 Decays into neutral Higgs bosons 

Now we turn to the decays involving neutral Higgs bosons. The channels Xa Xj^-fc (j = 
1,2,3; k = 1,2,3) can serve as sources for Higgs production from SUSY cascades at the 
LHC, and are therefore of particular interest. 

The decay Xa ~^ Xi^i is shown in Fig. [91 Contrary to what we observed for the decays 
into charginos, the two scenarios result in very different decay widths and both show a strong 
dependence on (pMi, with partial widths varying between 0.45 GeV and 0.13 GeV for Sh and 
between 0.06 GeV and 0.01 GeV for S g . The strong dependence on the phase <fMi is a 
consequence of the change in the relative CP-phase of Xa an d Xii while the corresponding 
CP-parity of the light Higgs, which is CP-even here, is not strongly affected. For ifM 1 = 0, 
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fMi ¥>Mi 

Figure 8: r(x° — >• XiW + ). Tree- level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [TJ with varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 



Xi and Xi have the same relative CP-parity, while for tfMx = ft they have the opposite one. 
Therefore, at (f>M 1 = 7r the decay is p-wave suppressed while at (f>Mi = the s-wave mode is 
allowed. In S h the partial decay widths are only partially suppressed, due to the relatively 
large phase space. In S g the suppression is stronger for the tree-level amplitude, leading 
to larger relative corrections. The relative corrections, shown in the upper right panel, are 
~ 10% for Sh and are between 2% and 8% for S g . For S g we observe a small difference 
between the two schemes at <fMi ~ ±90° , of the order of 0.1% (see also Tab. E]). It should 
be noted that here, as opposed to the case of the previous decays, the fact that the tree- 
level decay width depends strongly on <pMi leads to a noticeable difference between the two 
schemes. This difference has been highlighted in Fig. [271 In the lower left panel we show 
the branching ratios and in the lower right panel its relative corrections. Since the difference 
between the schemes is here very small we only show these results in scheme I. The peaks 
at (fMi = 35° and 325° are due to the threshold for the decay x® ~~ >" Xi^i, which leads to a 
singularity which affects the total width (see the discussion below on these threshold effects). 
It should be noted that the decay widths and the corresponding branching ratios, as well 
as their relative corrections, are roughly proportional here because the total width of x% 
shown in Fig. [2H below, is almost independent of <£mi m both scenarios, see the discussion 
in Sec. 14.61 

The results for the decay xi ~~ ^ Xi^2 are shown in Fig. [TUJ Since in both scenarios h,2 
tends to the CP-odd Higgs boson for real couplings, the dependence on <fMi is opposite to 
that for the decay x° — ► Xi^-i discussed above, with a p-wave suppression at <p Ml = and 
with the s-wave mode being allowed at 0m 1 = n. The decay widths are a factor two smaller 
than those for the decay into the lightest Higgs for Sh, and of the same order for S g . The 
relative corrections are also similar, of order ~ 10% in Sh and between 2% and 10% for S g . 
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Figure 9: ^(xl ~~ >* Xi^i)- Tree-level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [T], with varied. The 
upper left plot shows the decay width, the upper right plot shows the relative size of the 
corrections. The lower left plot shows the BR, the lower right plot shows the relative size of 
the BR. 



We can also observe a small difference at large ip Ml between the two schemes in S g , of the 
order of 0.05%, (see Tab. [5]). 

The decay xl — > X?^3 is shown in Fig. [TTJ The ip Ml dependence of the partial width and 
its relative correction for this process is qualitatively similar to that of xl ~~ >* Xi^i- However, 
for Sh the partial width is much smaller, between 0.01 GeV and 0.04 GeV. As the difference 
between the two schemes is of the order of 0.01%, (see Tab. EJ), they can not be visibly 
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Figure 10: r(x° — >• Xih 2 ). Tree-level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [TJ with (fM 1 varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 

T[GeV] AT/T[%] 




<PMi <PMi 



Figure 11: r(x° — > Xih 3 ). Tree-level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [H, with (pu\ varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 



distinguished here. 

The decays xt ~^ X^i, x® ~~ >* xH^i Xa ~~ >" X2^3 are shown in Figs. [TJl fT31and[TH respec- 
tively. In scenario Sh the second lightest neutralino is mainly wino-like, with a small mixing 
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Figure 12: r(x° — > x^hi). Tree-level ("tree") and full one- loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [TJ with (fM 1 varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 



with the bino component, leading to a weak dependence on tpu x - For S g the dependence 
on tpM 1 is much larger since the second lightest neutralino has both large Higgsino and bino 
components. 

For x°a ~^ X2^-i5 i n as it is shown in Fig. [121 the decay width oscillates from ~ 1.05 
to 1.2 GeV, and for Sh the decay width is ~ 1.2 GeV. It should be noted that this decay 
is s-wave mode allowed, and therefore the decay widths are 2.5 to 9 times larger than the 
decay into the lightest Higgs and neutralino. This turns out to be the dominating process, 
with branching ratios of up to ~ 8% and ~ 19% for, respectively, S g and Sh- The larger 
ipMr dependence in S g is due to the strong bino-Higgsino mixing of x\ in this scenario. This 
feature will be equally relevant for the remaining decays into either x\ or X3 discussed below. 
For S g , we observe the effect of the threshold for x\ — > Xi^i a t V^Afi = 35° and 325°. The 
dips are due to the resulting singular behavior of the derivatives of the self energies entering 
the field renormalization constants. This effect will also be observed in the other decays to 
a see Figs. [T3J, [TU and [191 i n the total width of Fig. [2U and in the branching ratios, 
see e.g. Fig. [91 It should be noted that both schemes have the same dips at <fMi = 35° and 
325°. This will be true as well for the other decays to X2 described below. The corrections 
are relatively small, ~ 1% to —3% for Sh and ~ 2% for Sh- In Fig. [T21 the renormalization 
schemes cannot be visibly distinguished from each other. The difference is below 0.01% for 
S g , and below 10~ 5 for Sh, see Tab. [5] for details. 

For X4 X2^2, in S g , as shown in Fig. [131 the decay width oscillates from ~ 0.03 to 
0.11 GeV, while for Sh it is ~ 0.03 GeV. The corresponding branching ratios in S g and Sh 
are, respectively, ~ 0.2 — 0.8% and ~ 0.5%. It should be noted that this decay is s-wave 
mode suppressed, and the decay width is therefore smaller than that for the decays into hi 
and /13. Again for S g the effect of the threshold for x\ ~^ at V^Mi = 35° and 325° is 
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Figure 13: r(x° — > x^)- Tree-level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [TJ with (fM 1 varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 



visible. The corrections are comparatively large for S g , ~ to 35%, and for Sh are ~ 6%. 
Despite a relatively significant difference between the schemes of 0.05% for S g , in Fig. CE3] 
this remains invisible. The large relative corrections in S g are a consequence of both the 
suppressed tree- level result, as well as the strong effect of the corrections on the mixing of 
the second lightest neutralino. 

The decay x°a X2^3> shown in Fig. [14j is s-wave mode allowed, and qualitatively similar 
to xl Xzhi- However, its decay width is smaller due to phase space suppression, oscillating 
from 0.14 to 0.18 GeV for S g and 0.12 GeV for Sh, and corrections are sizeable, ~ —5% and 
15% for S g and Sh, respectively. The two schemes cannot be visibly distinguished here. The 
difference between the schemes reaches 0.01% for S g , and again this remains invisible (see 
Tab. ED). 

The decays x°a - > X3h.i1 Xa ~^ xt^, X® ~^ xt^, shown in Figs. [T5lfT6land[T71 respectively, 
are kinematically closed in scenario Sh- For S g , there is a strong dependence on (pM 1 since X3 
has both large Higgsino and bino components. However, contrary to what we observed for 
the decays to X2 m <Sg; the decay to h 2 is s-wave mode allowed for cp Ml = 0, while the other 
two decays are suppressed. This is due to the opposite relative CP-parity of the x\ — xl pair 
relative to the X2 ~ Xa pair. 

For x\ ~^ xV 1 !-, as shown in Fig. [13, an oscillating behaviour similar to that for x° — > x®^ 
but even more enhanced (going from ~ 0.01 to 0.13 GeV) results in even larger relative 
corrections. In fact the suppression of the tree level is now larger, due to the smaller phase 
space for this decay, and mixing with the unsuppressed states of the neutralinos has a 
dramatic effect, with corrections approaching 70%. Notice, however, that the second and 
third lightest neutralinos are roughly degenerate in S g , see Tab. HJ This leads naturally to a 
large mixing character for these two mass eigenstates, supported by the fact that the sum of 
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Figure 14: r(x° — > X2/13). Tree-level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [TJ with (fM 1 varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 



the decay widths is much less sensitive to (pMi- Therefore, the large corrections should not be 
regarded as a breakdown of the renormalization procedure but rather as an indication that 
one should consider all the neutralino states simultaneously. For this particular decay the 
branching ratio does not reach 1%. In Fig. the difference between the renormalization 
schemes reaches ~ 0.02% and they cannot be visibly distinguished from each other, see 
Tab. [5] for details. 

For the decays into the heavier Higgs bosons the corrections are mild, at the level of a few 
percent. In x°a — > xt^, as shown in Fig. (THl the unsuppressed decay, with widths between 
~ 0.6 and 0.7 GeV, receives small corrections via mixing with the p-wave suppressed states. 

In xl X3^3; as shown in Fig. [TTJ, the <j>M\ dependence is small due to the combination 
of one p-wave suppressed amplitude with an s-wave allowed one in which the couplings are 
small, resulting in corrections of a few percent. For these decays the small difference between 
the two schemes cannot be observed in the figures. 

4.3 Decays into Z bosons 

The channels involving the Z boson, x\ ~~ >* X^jZ-i are presented in Figs. [TSUiZUl The strong 
resemblance in the (fMi~ dependence between these plots and those for the decay into h 2 is 
due to fact that gauge bosons are CP-odd, while in our scenarios the second Higgs boson has 
a very small CP-even component and tends to the CP-odd state for ip^i — 0, vr. There is a 
strong dependence on the relative CP-phase of the neutralinos in the initial and final states, 
which lead to visible differences between the renormalization schemes, given explicitly in 
Tab. El 

In both scenarios the lightest neutralino is mainly bino-like, therefore xl ~~ ^ Xi% depends 
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Figure 15: r(x° — > x^/ii)- Tree-level ("tree") and full one- loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [TJ with (fM 1 varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 
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Figure 16: r(x° — > xjj/^)- Tree-level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [H, with (p^i varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 



strongly on <p Ml , as it is shown in Fig. [THJ The decay is both qualitatively and quantitatively 
similar to the case of x° — > X?^2 ; due to the relative CP-phase of the neutralinos, and the fact 
that the Z is CP-odd, the decay is suppressed at <p Ml = 0, and maximal at cp^ = 180°, and 
the decay width ranging from to 0.14 GeV and 0.08 to 0.25 GeV for S g and Sh respectively. 
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Figure 17: r(x° — > X3/13). Tree-level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [TJ with (fM 1 varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 




Figure 18: T(x° — > X\Z)- Tree-level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [H, with y>M\ varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 



For S g the dependence of the relative size of the corrections on (fM x is much larger, from 1% 
to 15%, with visible differences between the schemes at large <£>Mi of up to 0.2%, while the 
loop corrections for Sh are almost independent of (fMu of ~ 8%. The difference between the 
schemes in S g has been highlighted in Fig. [271 
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Figure 19: r(x° — > x^)- Tree- level ("tree") and full one- loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [TJ with ^>m x varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 



For x\ in S g the decay width oscillates from ~ 0.05 to 0.15 GeV and for Sh 

the decay width is ~ 0.06 GeV, as it is shown in Fig. [19j For S g the effect of the threshold 
for X2 ~* Xih-i a ^ <PMi = 35° and 325° is visible as a small dip in the decay width and a 
marked dip in the relative corrections. The corrections are comparatively large for S g , ~ 
to 25%, and for Sh are ~ 3%. The two schemes are not visibly distinguished from each 
other, with the largest differences in S g of ~ 4 x 10~ 4 . Although in Sh it is below threshold, 
in S g the decay x® ~^ X3Z, shown in Fig. [2Q1 is s-wave mode allowed in S g , resulting in the 
largest branching ratio into Z, while it is below threshold in Sh- The decay width goes from 
0.78 GeV at <fMi ~ to 0.62 GeV, with the corrections between ~ 6.2 and 10.5%. The two 
schemes differ by up to ~ 5 x 10~ 5 and cannot be visibly distinguished. 

4.4 Decays into (s)leptons 

Now we turn to the decays involving (scalar) leptons. The expressions for all these decay 
widths follow the same pattern, see the expressions for the tree level widths in Appendix [A] 
The dependence on ip^ is small, although the results for Sh do show some dependence due 
to the small bino-like component of the decaying neutralino. We have chosen M~ 1l < Mg R , 
leading to lighter left-handed and heavier right-handed sleptons, and significant mixing in 
the scalar tau sector. 

In Fig. [21] we show the results for the decay x°a ~* T ~ ^1 ■ The decay widths are found 
to be an order of magnitude larger in S g (~ 0.37 GeV) than in Sh (~ 0.03 GeV), since the 
gaugino-like neutralino has an unsuppressed coupling to the large left component of the stau, 
while the Higgsino-like neutralino couples to the suppressed Yukawa coupling. This pattern 
is even more significant for the decays into the lower generation sleptons not shown here. In 
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Figure 20: r(x° — > X3Z). Tree- level ("tree") and full one- loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [TJ with ^>m x varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 




Figure 21: r(x° — > t + t{). Tree-level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [TJ with varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 



the right panel we observe that the one- loop corrections are very small in S g , while they are 
around ~ 20% in Sh, due to the larger dependence on the stau mixing 

The results for the decay into the heavier scalar tau are shown in Fig. |22j The pattern 
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Figure 22: r(x° — >■ r + f 2 ). Tree-level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. [T], with varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 



is similar to the preceding decay, the difference being that here the right-handed component 
of the heavier stau is larger, resulting in a decay width which is four times larger for Sh, and 
30% smaller for S g . The corrections in S g remain very small while those in Sh are now ~ 7%. 
For the decay into the first two generations, where the slepton mixing is usually negligible, 
both scenarios have very small partial widths. 

The results for the decay x\ — ^ v t are shown in Fig. [221 Here, the gaugino-like 
neutralino has a decay width which is roughly the sum of decay widths for x° — > t~ and 
xl —> r~ f£, i-e. ~ 0.7 GeV for S g and ~ 0.05 GeV for Sh- The radiative corrections are 
~ 0.5% in S g and ~ —6% in Sh- 
in S g , where the neutralino is gaugino-like, the branching ratios for these tree processes 
are, respectively ~ 2.5%, ~ 2%, and ~ 5%. Taking into account the charged conjugated 
processes this results in a branching ratio of almost 20% for the third lepton family. For 
the first two generations the branching ratios into every left-handed slepton or sneutrino is 
~ 4.5 and 5%, respectively, while the decays into the right-handed sleptons are negligible. 
Therefore, in this class of scenarios the leptonic decays could be the dominant ones. In Sh, on 
the other hand, the leptonic decays are subdominant, especially for the first to generations, 
and mainly due to the small gaugino component of the decaying neutralino. 

In all the decays to leptons, the difference between the renormalization schemes is negli- 
gible, as shown in Figs. [211 [22] and [231 an d summarized in Tab. [5j 

4.5 Full one-loop results: total decay widths 

In this subsection we briefly show the results for the total decay widths and its relative 
correction of the three heaviest neutralinos in S g and Sh- The decay width of x% shown in 
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Figure 23: r(x° — > u T u\). Tree-level ("tree") and full one-loop ("full") corrected decay 
widths are shown with the parameters chosen according to Tab. dj with varied. The 
left plot shows the decay width, the right plot shows the relative size of the corrections. 



the l.h.s. of Fig. [22J is almost independent of <pMi, as is expected from the heavier neutralino 
in a GUT-related scenario. Therefore, the branching ratios can be easily obtained from the 
partial widths. The large (pMx dependence of the single channels is due to the strong effect 
on the mixing of the different neutralinos, as already argued in the preceding subsections. 
The corrections are also very small, of ~ ±1%, with the expected dips due to the thresholds 
for x°2 -> Xih at <p Ml = 35° and 325° and Xi ->■ Xx^~ at <p Ml = 125° and 235°. 

The total decay width of x® is shown in the l.h.s. of Fig. |25j As already discussed, in 
S g the second and third lightest neutralinos have similar masses and have a Higgsino-bino 
mixing character which strongly depends on pu x - This is both true at the tree-level, where 
the total width goes from ~ 0.04 GeV for <p Ml =0 to almost three times as much for ip^ = 7r, 
as well as for the loop corrections, which are of O(20%). The total width is also significantly 
smaller than that of x1> largely due to the reduced phase space, see Tab. [2j Here only the 
leptonic decays and those to Z and the lightest neutralino are open. On the contrary, in 
Sh both X3 an d X°a are Higgsino-dominated and nearly degenerate. Consequently the same 
decay channels are open and their widths are similar, both at tree-level and at one-loop. 

The total decay width of xH ls shown in the l.h.s. of Fig. [26j The strong mixing of the 
second and third lightest neutralinos in S g has already been discussed. The same decay 
channels as for X3 are open, except for that to Z and the lightest neutralino which is only 
open for \<Pmi \ < 32°. However, the decay to Z is subdominant with a BR smaller than 1%. 
The threshold effect we observe in the right panel for the relative corrections is due to the 
decay into the lightest Higgs boson, as already discussed in this section for those decays with 
final x_2- The decay width and its relative correction show a complementary behavior to 

7 It should be noted that this effect, due to the singularity of the wave function renormalization, is 
characteristic of on-shell renormalization schemes, which are less precise when thresholds of external particles 
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Figure 24: F tot (x1). Tree- level ("tree") and full one-loop ("full") corrected total decay 
widths are shown with the parameters chosen according to Tab. [TJ with ^>m x varied. The 
left panel shows the decay width, the right panel shows the corresponding relative size of 
the corrections. 
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Figure 25: r tot (xg). Tree- level ("tree") and full one-loop ("full") corrected total decay 
widths are shown with the parameters chosen according to Tab. [Tj with varied. The 
left panel shows the decay width, the right panel shows the corresponding relative size of 
the corrections. 



are open. The masses entering these thresholds are the tree-level ones, which in the case of the lightest Higgs 
boson is close to that of the Z boson. The renormalized lightest Higgs boson, on the other hand, has a mass 
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Figure 26: ^tot(x%)- Tree- level ("tree") and full one-loop ("full") corrected total decay 
widths are shown with the parameters chosen according to Tab. [TJ with ^>m x varied. The 
left panel shows the decay width, the right panel shows the corresponding relative size of 
the corrections. 



the corresponding ones for the third lightest neutralino, i.e., the dependence on ^>m x of the 
sum of both widths is much weaker, with corrections of ~ 10%. 

In Sh, again, the neutralinos do not strongly mix and x® * s wino-dominated. Here the 
decays into left-handed sleptons or sneutrinos dominate due to the strong wino coupling. 
However, the subdominant decays to the lightest neutralino and Z or Higgs bosons become 
the only open channels if the sleptons are chosen heavier. These decays show a strong (and 
complementary) dependence on ip Ml due to the change in the relative CP-parity of the two 
lightest neutralinos. 

4.6 Differences between the renormalization schemes 

In our benchmark scenarios S g and Sh we have found remarkably good agreement between 
the two schemes. For most decay channels, the difference between the relative corrections to 
the partial widths AT/T for scheme I and II is of the order of 10~ 5 . In Tab. |5]we show this 
difference for all the decay channels, in both S g and Sh, and find that the largest differences 
are observed in xl ~> xl -> xl h 2, xl -> xih, xl -> Xi,2 Z - In Fig. [27] we highlight 

these differences in S g for xl ~~ >* Xi z and x° — * Xi^x- To be precise, in this figure we 
compare the full one-loop correction Ar calculated in scheme I both with and without the 
squark shifts (see end of Sec. 12. 21) to that calculated in scheme II where squark shifts are not 
included. One can see that at — 0° and 180° the difference between the results without 
squark shifts vanishes, confirming that the schemes differ only in the treatment of the phases. 
One can also clearly see the impact of the squark shifts on the size of the one-loop correction. 

of ~ 126 GeV and the decay is closed. 
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Channel 



xl -»> XiH + 
xl XiW + 

xl -»> x?/n 

X4 -> X?^2 

X° ->■ Xih 3 
Xl X^i 
X^ -»> xjjto 



X4 -> X^S 

X° -> Xs^i 
X° -»• X° 3 h 2 
Xl X^3 
X^ -X X?^ 
X^ ^ X^ 
X^ -> X^ 

X4 -> T- f+ 

Xl r" f + 

X4 ~> g 



45° 



1.5 x 1(T 5 

3.4 x 1(T 6 
-1.9 x 10~ 4 

4.5 x 10~ 4 
-1.4 x 10~ 4 
-1.2 x 10~ 5 

1.3 x 1(T 4 
5.5 x 10~ 5 

2.4 x 10~ 4 
-2.8 x 1(T 5 
-5.1 x 10" 5 

1.2 x 1(T 3 
2.2 x 10" 4 
-4.5 x 10" 5 
-1.3 x 1(T 6 
-6.5 x 10" 5 
-3.5 x 1(T 6 



90° 



2.0 x 10~ 5 

6.1 x 10~ 6 
-6.1 x HT 4 

5.2 x 10~ 4 
-3.6 x 10~ 4 

5.4 x HT 5 
-4.9 x 10~ 4 

1.0 x 10~ 4 
-7.2 x lO" 4 

3.2 x 10~ 5 
-1.6 x 10~ 4 

7.5 x 10^ 4 
-3.4 x HT 4 

1.2 x 10~ 5 
3.7 x 10~ 6 
-6.2 x 10~ 5 
-7.0 x 10~ 6 



S h 



45 c 



-6.9 x 10~ 6 
9.9 x 10~ 5 

-6.3 x 10" 5 
1.5 x 10~ 4 

-9.1 x 10~ 5 
2.7 x 10~ 6 
6.0 x 10~ 6 
4.5 x 10- 6 



1.2 x HT 4 
9.5 x 10~ 6 

-3.5 x 10" 5 
-1.1 x 10~ 4 
-8.9 x 10~ 6 



90° 



-4.8 x 10~ 6 
9.7 x 10~ 5 

-1.8 x 10~ 4 
1.7 x 10- 4 

-2.2 x 10~ 4 
6.1 x 10~ 6 
7.1 x 10~ 6 
8.6 x 10- 6 



1.7 x 10~ 4 
7.3 x 10~ 6 

-1.1 x lO" 5 
-1.1 x lO" 4 
-1.9 x 10" 5 



Table 5: Differences between the relative corrections to the decay width for scheme I and 
II, shown in both scenarios S g and S h , i.e. AT/T(S g /S h ) — Ar/r(«S fl /<S^), at the specified 
values of (p Ml - The missing results correspond to those channels for which the decays are 
below threshold S h . 



Earlier in this section, it was noted that larger differences between the schemes are observed 
when the tree- level decay width depends strongly on the phase </?Mi- This explains why 
decays to the lightest neutralino are the most strongly affected ones. Also, the differences 
are in general much more pronounced in S g than in Sh, due to the mixing between the 
bino and Higgsino components, except for the decays into sleptons or a W^-boson, where the 
Higgsino component has a negligible role. 
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Figure 27: T(x4 — > Xi^) an d r(x° ~~ Xi^i)- The one-loop ("full") correction to the partial 
decay widths are shown as a function of (fMi- The parameters are chosen according to Tab. [TJ 
Also shown is the correction without including the shifts in the squark sector ("full no A<f' ). 



5 Conclusions 

We have evaluated all non-hadronic two-body decay widths of neutralinos in the Minimal 
Supersymmetric Standard Model with complex parameters (cMSSM). Assuming heavy scalar 
quarks we take into account all decay channels involving charginos, neutralinos, (scalar) 
leptons, Higgs bosons and SM gauge bosons. The decay modes are given in Eqs. (jSJ) - ©• 
The evaluation of the decay widths is based on a full one-loop calculation including hard 
and soft QED radiation. Such a calculation is necessary to derive a reliable prediction of any 
two-body branching ratio. Three-body decay modes can become sizable only if all the two- 
body channels are kinematically (nearly) closed and have thus been neglected throughout 
the paper. The same applies to two-body decay modes that appear only at the one-loop 
level. 

We first reviewed the one-loop renormalization of the cMSSM, concentrating on the most 
relevant aspects for our calculation, except for the details for the Higgs boson sector which 
can be found in Ref. [16]. More importantly, we have given details for the chargino/neu- 
tralino sector in the two on-shell renormalization schemes which we have compared in this 
work. The two schemes differ in the treatment of complex contributions in the chargino and 
neutralino sector. The different renormalization of the CP-violating phases leads to small 
differences in the cMSSM, which are however of higher order in the electroweak coupling and 
vanish in the limit of real couplings. Differences indicate the size of unknown higher-order 
corrections involving complex phases beyond the one-loop level. We have also discussed the 
calculation of the one-loop diagrams, and the treatment of UV- and IR-divergences that 
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are canceled by the inclusion of soft QED radiation. Our calculation set-up can easily be 
extended to other two-body decays involving (scalar) quarks. We have taken into account 
all absorptive contributions, explicitly inlcuding those of self-energy type on external legs. 
This ensures that all CP-violating effects are correctly accounted for. 

In the numerical analysis we mainly concentrated on the decays of the heaviest neutralino, 
xl- For this analysis we have chosen a parameter set that allows simultaneously all two-body 
decay modes under investigation, and respects the current experimental bounds on Higgs 
boson and SUSY searches (where the combination of low M H ± and relatively large tan (3 is in 
a potential conflict with the most recent LHC searches for the heavy MSSM Higgs bosons). 
The masses of the charginos, and thus roughly those of the second and fourth neutralino, are 
in this scenario, respectively, 350 and 600 GeV. This leads to two representative scenarios 
for the chargino/neutralino sector, for fi > M 2 or \i < M 2 . These benchmark scenarios 
allow copious production of the neutralinos in SUSY cascades at the LHC. Furthermore, 
the production of xlx°j at the ILC(IOOO), i.e. with y/s = 1000 GeV, via e+e" ->• xlx°j will 
be possible, with all the subsequent decay modes (J3j) (J6]) being (in principle) open. The 
clean environment of the ILC would then permit a detailed, statistically dominated study 
of the neutralino decays. Depending on the channel and the polarization, a precision at the 
percent level seems to be achievable. Special attention is paid to neutralino decays involving 
the Lightest Supersymmetric Particle (LSP), i.e. the lightest neutralino, or a neutral or 
charged Higgs boson. 

We have shown results for varying <£>Mi, the phase of the soft SUSY-breaking parame- 
ter Mi, which leads to CP- violation in the chargino and neutralino sectors. We have analyzed 
the tree-level and full one-loop results for all kinematically open decay channels of the heav- 
iest neutralino. For the decays of the second and third neutralinos we have only shown the 
total widths. 

We found sizable corrections in many of the decay channels. The higher-order corrections 
of the neutralino decay widths involving the LSP are generically up to a level of about 10%, 
and decay modes involving Higgs bosons can easily have corrections up to 20 — 30%. The size 
of the full one-loop corrections to the decay widths and the branching ratios also depends 
strongly on (p^n especially for those decays in which an external neutralino is a mixed 
Higgsino-bino state. We conclude that the largest effect of the radiative corrections is due 
to its effect on the mixing of the neutralinos. All results on partial decay widths of x° as 
well as the total decay widths of all neutralinos are given in detail in Sec. HI 

For the two on-shell renormalization schemes considered, we have found very good agree- 
ment: where the difference between the relative size of the corrections is found to be < 0.1%. 
The largestest differences have been found for decays with a strong dependence on the pa- 
rameter (fMi- The good agreement between the two schemes is not unexpected, as they are 
found to be equivalent up to a higher order effect, as discussed in Sec. 14.61 

The numerical results we have shown are of course dependent on choice of the MSSM 
parameters. Nevertheless, they give an idea of the relevance of the full one-loop corrections. 
For other choices of SUSY masses the corrections to the decay widths would stay the same, 
but the branching ratios would look very different. Channels for which the decay width 
(and its respective one-loop corrections) may look unobservable due to the smallness of the 
BR in our numerical examples could become important if other channels are kinematically 
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forbidden. 

Following our analysis it is evident that the full one-loop corrections are mandatory for 
a precise prediction of the various branching ratios. This applies to LHC analyses, but even 
more to analyses at the ILC or CLIC, where a precision at the percent level is anticipated for 
the determination of neutralino branching ratios (depending on the neutralino masses, the 
center-of-mass energy and the integrated luminosity). The results for the neutralino decays 
will be implemented into the Fortran code FeynHiggs. 
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Appendix 



A Tree-level results 

For completeness we include here the expressions for the tree-level decay widths: 

r tree (x° -X XjH+) = [(\C(xlxl,H-) L \ 2 + \C(xlxl,H-) R \ 2 ) (m\,+m^-Ml ± ] 
+4 Re {C{xl Xt H-)lC(xl X- , H~) R } m^i 



,m xt 



x 



327rm 3 



(2 = 2,3,4, j = 1,(2)) , 



r tree (x° x~w + ) = [(\c(xl xf, w~) L \ 2 + \c(xl x}> w~) R \ 2 ) 



(ml — m 2 ± ) 2 
x|4 + ^-2M^^^- 



-12 Re {C(x°, xj, W-)* L C(xl X+ W~) R } m^om-± 



x 



327rm 3 



(2 = 2,3,4, j = 1,(2)) , 



T tree (x° -X X]h k ) = [{\C{xl X% h) L \ 2 + \C(xl X% h k ) R \ 2 ) (m| ? + 
+4 Re {C(xl X% h k )* L C(xl X°, h k ) R ) 



Xj X j 



m l k ) 



(A.I) 



(A.2) 
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X 1 / 2 (m 2 ,m 2 ,m 2 hk ) 



x 



32nm% 

Xi 



2,3,4, j = 1,(2), k = 1,2, 3) 



+2 Re 



x 



A 1/2 (m 2 ,m 2 ,m 2 



167rm_ 



(z = 2,3,4, j = 1,(2), A; = 1,2,3) 



-itree /• ~0 



(m 2 - m 2 ) 2 " 
x I m 2 + m~o - 2M 2 + - x 



A' 



Ml 



Xi X 7 



A 1/2 (m| , m| ,M§) 

X ^ 

327rm , i 



[|C(x°,x;,^) R | 2 m 2 o+m 2 -2M 2 + 



+6 Re 



(m 2 -m 2 ) 2, 

Ai A.j 



x 



A 1/2 (m 2 o,m 2 ,M 

A,i Xj 



X Q i,XpZ) R ) Jm^m^o 

(z = 2,3,4, j < i) 



167rm? 



r tree ( ~U _^ £ - ^ = ^ ^ ^ + | C( ~U ? £? ( , + ^ _ , j 



+4 Re [C(x°, I ik)lC(xl I Qr\ m%m t 
A 1 / 2 (m|o,m^,m|^ 



r tree (x° -> ^ ^ T ) = \C(xl n, n) R \ 2 {mL - m 2 ) 



327rm? 



i = 2,3,4, 



e,/z,r, fc = 1,2) 



x 



A 1/2 (m? ,0, m| 
Xi 



327rm? 



2,3,4, £ = e,/i,r) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



where A(x, y,z) = [x — y — z) 2 — 4yz and the couplings C(a,b,c) can be found in the 
FeynArts model files jH2]. C(a, b, c) l.r denote the part of the coupling which is proportional 

to uj t = |(i T75 ). 
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